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Abstract 

The replication of many RNA viruses involves the translation of polyproteins, whose processing by endopeptidases is a 
critical step for the release of functional subunits. PI is the first protease encoded in plant potyvirus genomes; once 
activated by an as-yet-unknown host factor, it acts in cis on its own C-terminal end, hydrolyzing the PI-HCPro junction. 
Earlier research suggests that PI cooperates with HCPro to inhibit host RNA silencing defenses. Using Plum pox virus as a 
model, we show that although PI does not have a major direct role in RNA silencing suppression, it can indeed modulate 
HCPro function by its self-cleavage activity. To study PI protease regulation, we used bioinformatic analysis and in vitro 
activity experiments to map the core C-terminal catalytic domain. We present evidence that the hypervariable region that 
precedes the protease domain is predicted as intrinsically disordered, and that it behaves as a negative regulator of PI 
proteolytic activity in in vitro cleavage assays. In viral infections, removal of the PI protease antagonistic regulator is 
associated with greater symptom severity, induction of salicylate-dependent pathogenesis-related proteins, and reduced 
viral loads. We suggest that fine modulation of a viral protease activity has evolved to keep viral amplification below host- 
detrimental levels, and thus to maintain higher long-term replicative capacity. 
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Introduction 

Viruses are obligate parasite pathogens that hijack host factors 
to assure their own survival and propagation. Due to the limited 
coding capacity of their genome, viruses undertake distinct 
translational strategies [1]; one of most widely employed by 
RNA viruses involves polyprotein synthesis [2]. A full-length 
polyprotein is theoretically derived from a single translational 
event, compromising the timely expression of the individual viral 
cistrons. To overcome this possible drawback and successfully 
regulate replication, assembly and spreading stages, various post- 
translational mechanisms have evolved to modulate the spatial- 
temporal availability of functional viral subunits. For instance, it 
is not uncommon that the same polyprotein is hydrolyzed by 
several endopeptidases; cleavage kinetics are thus linked to 
enzyme processivity and, in trans-acting proteases, to the different 
alTmity with the specific cleavage sites [3,4]. Activation of viral 
proteases might depend on the availability of defined cell- or 
pathogen-encoded cofactors [5-7] and structural rearrangements 
that modulate substrate accessibility, as shown for the hepatitis C 
virus NS3 protease domain [8]. Zymogen activation and 
allostery are key regulatory mechanisms of trypsin-like proteases 
[9], a group of enzymes that is widespread in positive strand 
RNA viruses [10] and that includes PI proteins of potyviruses 
[11]. 



Members of the genus Potyvirus (family Potyviridae) belong to the 
picorna-like supergroup and represent one of the largest groups of 
plant-infecting RNA viruses [12,13]. Their single-stranded RNA 
genome is ~ 1 0 kb in size and encodes a large polyprotein 
comprising (from N- to C-terminus) PI, HCPro, P3, 6K1, CI, 
6K2, VPg, NIa-Pro, Nib, and the coat protein (CP). An additional 
protein, P3N-PIPO, is originated by a frameshift in the P3 cistron 
[14,15]. The cysteine protease NIa-Pro processes the C-terminal 
part of the polyprotein by seven cleavage events, while PI and 
HCPro are responsible for their own release by a cis cleavage at 
their respective C-termini [16-18]. Once released, the mature PI 
and HCPro C-terminal extremities are thought to be trapped in 
the active cleft, leading to autoinhibition of trans cleavage activity 
[18,19]. Located at the beginning of the polyprotein, PI was the 
last potyviral endopeptidase identified; inactivating mutations of its 
catalytic domain preclude virus viability [20], making PI an 
attractive target for the development of antiviral tools. In contrast 
to the two other genome-encoded proteases, PI relies on a still 
unidentified host factor for its activation [18]. Computational 
analysis of PI potyviral proteins showed its great variability both in 
length and in amino acid sequence, and its diversification in 
potyviral species was thus associated with host specialization [21]. 
Although PI involvement in the definition of virus host range 
was higlrlighted [22,23], its specific contribution to potyviral 
infection is still unclear. Many functions were attributed to P 1 , such 
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Author Summary 

RNA viruses are ideal systems for tlie study of population 
dynamics, relationships among pathogen traits such as 
fitness and virulence, and of host immune responses to 
pathogen attacks. Based on experimental evolution stud- 
ies, early models equated parasite virulence with fitness. 
Some reports showed that viral virulence and fitness can 
be unlinked. Here we present evidence that the highly 
disordered N-terminal region of a potyviral PI protein 
negatively regulates its self-cleavage activity. Removal of 
this regulator domain greatly affects viral infection, which 
is characterized by accelerated early replication and 
enhanced symptom severity. These properties are none- 
theless associated with low viral accumulation and high 
induction of antiviral resistance markers. Finally, we 
propose that host-dependent regulation of PI processing 
efficiency modulates viral virulence and alleviates the host 
antiviral responses. 
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as cell-to-cell movement, systemic spread, and viral genome 
replication enhancement [reviewed in 24]; PI was later shown to 
strengthen the RNA silencing suppressor activity of HCPro [25-28]. 

Here we used Plum pox virus (PPV) as a model system to study 
potyviral PI. Our results indicate that P 1 and its protease domain 
have a positive role in potyvirus infection, independent of silencing 
suppression. Moreover, the in vitro and inplanta data presented here 
suggest that host-dependent regulation of PI self-cleavage activity 
has evolved to modulate the efficiency of viral infection and escape 
plant defense responses. 

Results 

PI deficiency in PPV Infection is not complemented by 
Arabidopsis plants with defective RNA silencing pathways 

Genetic rescue of a Turnip mosaic potyvirus (TuMV) defective in 
RNA silencing suppression activity was reported in the 
Arabidopsis thaliana triple mutant line dcl2-l dcl3-l dcl4-2 {dd2/ 
3/4) [29]. A PPV isolate adapted to Micotiana develandii and able 
to infect herbaceous hosts [30], as well as its infectious cDNA 
clone encoding a green fluorescent protein (GFP) reporter gene, 
were described [31]. To study the PI contribution in silencing 
suppression, we generated a viral clone that lacks the entire PI 
sequence (API, Figure lA), which we used to challenge 
Arabidopsis mutant plants with defective antiviral silencing 
pathways. At 15 days post-agro-inoculation (dpi), API -infected 
Col-0 and dcl2/3/4 plants were almost symptomless, whereas 
PPV-infected Col-0 plants showed leaf chlorosis, which was 
stronger in the dcl2l3 / 4 mutant line (Figure IB). The presence of 
API and wild-type PPV in systemic leaves was evidenced by GFP 
fluorescence detection (Figure IB). Anti-PPV CP western blot 
analysis of systemicaUy infected leaves showed that, in both host 
genotypes, viral accumulation was significantly higher in wild- 
type PPV than in API -infected plants (Figure IC and D). The 
results confirmed that lack of PI sequence did not compromise 
virus ability to replicate or to move systemically, as also shown 
for another potyvirus, Tobacco etch vims (TEV) [32]. Although 
API appears to replicate better in the silencing deficient plants 
than in wild-type Col-0, the dcl2/3/4 line did not restore the 
API phenotype and viral accumulation to wild-type PPV 
levels (Figure IB and D). Based on these results, the main role 
of P 1 in viral amplification appears unrelated to RNA silencing 
suppression. 
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Figure 1. Arabidopsis mutant plants with defective RNA 
silencing pathways fail to rescue PPV AP1 amplification 
defects. (A) Representation of wild-type PPV and API, a PPV clone 
with deletion of the entire PI sequence. The reporter sGFP(S65T) [103] 
gene is present between Nib and CP coding sequences of both viral 
clones. Boxes with diagonal lines indicate P3N-PIP0 protein. (B) 
Symptoms and GFP fluorescence of Arabidopsis Col-0 and dcl2/3/4 
plants agro-inoculated with wild-type PPV or with APT Pictures were 
taken in an epifluorescence microscope at 15 dpi. Scale bar, 1 cm. (C) 
API and wild-type PPV viral accumulation in A. thaliana Col-0 systemic 
leaves was assessed by anti-PPV CP (CP) western blot assay (15 dpi). 
Each lane represents a pool of samples used for signal quantification. H, 
healthy Col-0 plant sample. Relative CP signal intensities were plotted 
using average PPV value equal to 100. (D) API and wild-type PPV viral 
accumulation in A. thaliana dcl2/3/4 systemic leaves was assessed by 
anti-PPV CP (CP) western blot assay of upper non-inoculated leaves 
(15 dpi). Each lane represents a pool of samples used for signal 
quantification. H, healthy dcl2/3/4 plant sample. Relative CP signal 
intensities were plotted using average PPV value equal to 100. 
Histograms show mean ±SD (n = 6 samples/condition, from two 
independent Agrobacterium cultures); *** p<0.001. Student's f-test. 
Ponceau red-stained blots are shown as loading control. 
doi:1 0.1 371 /journal.ppat.1 003985.g001 

The N-terminal region of potyviral PI is predicted to be 
intrinsically disordered 

To further analyze possible PI functions, we coupled mul- 
tiple sequence alignment and structural predictors for bioinfor- 
matic analysis of its sequence. Intrinsic protein disorder was 
estimated using DISOPRED2 [33] and MetaDisorderMD2 
[34]. 

From the in silico analysis of P 1 proteins from different potyviral 
species groups (Table SI, in Text SI), we identified FxxLE as a 
conserved motif in the PI N-terminal region (Figure 2), in 
conjunction with the reported IxFG and ISI motifs [21]. The 
relatively well-conserved C-terminal regions of the proteins, which 
correspond to the protease domain (defined below) and include the 
VELI motif, are predicted mainly as structured. The residues from 
N-terminal sequence patches with least conservation are generally 
predicted to be intrinsically unstructured (Figure 2), suggesting that 
final protein conformation is the main evolutionary target, rather 
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Figure 2. Alignment of PI amino acid sequences from PPV and six reference potyviruses. GenBank accession numbers are reported in 
Table SI, in Text SI; amino acid background is assigned according to the ClustalX color scheme [120], Residues aligning to PPV PI minimal protease 
domain (aa 1 65-308, based on our analysis) are shown in black letters and bright-colored background. Alignment quality, based on BLOSUIV162 scores, 
is shown as a bar graph (below). FxxLE and the previously reported conserved motifs IxFG, ISI, VELI are boxed. The protease conserved RG dipeptide is 
marked with black bar and the catalytic triad His, Asp and Ser is marked with stars. PPV IVID2 and PPV DP2 lines show protein disorder prediction of 
PPV PI sequence according to IVIetadisorderMD2 and DIS0PRED2 analyses, respectively, where "+" are disordered and "-" ordered residues. 
Additionally, D1S0PRED2 prediction confidence for each PI sequence is plotted: PPV, black line; PVY, orange; PVA, turquoise; PVV, magenta; TuMV, 
blue; WY, purple; YMV, green. The first 46 residues of TuMV PI are hidden, as they do not align to the other accessions. 
doi:10.1371/journal.ppat.1003985.g002 



than maintenance of the primary sequence, as proposed [35]. The 
disorder predictions were furtlier supported by PI sequence 
analysis of Scallion mosaic virus and the monocotyledon-infecting 



Sugarcane mosaic virus group (Figure SI), which comprise the smallest 
known PI sequences and are suggested to be at the base of the 
potyviral evolutionary tree [13]. 
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Figure 3. Involvement of P1 conserved motifs in protease 
cleavage. (A) Diagram of the DNA constructs used in the transcription 
reactions. T7 RNA polymerase promoter, black arrow, drives the cDNA 
of PPV 1-1361 genome nucleotides, comprising the 5'UTR and the ORF 
that codes for PI (black box) and 97 N-terminal amino acids from HCPro 
(HC-97, grey box). HC-97 is followed by a stop codon and PPV 3'UTR 
cDNA. PI residues chosen for alanine scanning are labeled (details 
below). When indicated, Thosea asigna virus 2A 'self-cleaving' peptide 
(T2A) was inserted in between PI and HC-97. Below, expected 
translation products and their molecular weights are shown. (B) PI 
self-cleavage was evaluated by in vitro translation in the wheat germ 
extract system. Labeled translation products were resolved by tricine- 
SDS-PAGE and the ^^S signal detected. WT, wild-type PI; FG, FG-6,7-AA 
substitution of the IxFG motif; WE, W63A E67A substitution of the FxxLE 
motif; VE, VE-189,190-AA substitution of the VELI motif; S, replacement 
by alanine of the catalytic S259; HY, HY-307,308-AA substitution of Pj 
and P, cleavage site residues; ST2A, catalytic S259A4-T2A. Molecular 
weight markers are indicated (left). 
doi:10.1371/journal.ppat.1003985.g003 

Involvement of PI conserved motifs in protease self- 
cleavage 

To study the relevance of PI motifs in protease activity, we 
mutated selected PPV PI conserved amino acids to alanine 
(Figure 3A). The viral cDNA constructs were transcribed and in 
ra'fro-translated using the wheat germ extract (WGE) system. 
Replacement of catalytic S259 (S) and P2 and Pi cleavage site 
residues HY-307,308 (HY with alanine impaired PI protease self- 
cleavage, precluding release of the mature 35.3 kDa PI protein 
and the 11.1 kDa HCPro fragment (HC-97) (Figure 3B). No 
proteolytic processing was detected in the construct bearing the 
VE-189,190-AA (VE) substitution in the VELI motif, indicating 
that it indeed belongs to the minimal protease domain, consistent 
with the disorder predictions. We found no appreciable 
differences between the wild-type PI (WT) construct and those 




Figure 4. GFP fluorescence and viral accumulation of PPV 
cDNA clones with altered PI sequence. N. ckvelandii plants were 
agro-inoculated with the indicated binary vectors: PPV, pSN-PPV; API, 
pSN-PPV API; FG, pSN-PPV P1-FG; WE, pSN-PPV PI -WE; VE, pSN-PPV P1- 
VE; S, pSN-PPV P1-S; HY, pSN-PPV P1-HY. (A) Pictures of upper non- 
inoculated leaves were taken in an epifluorescence microscope (22 dpi). 
Scale bar, 1 cm. (B) Anti-PPV CP (CP) western blot assay of agro- 
inoculated N. clevelandii leaves (10 dpi); each lane corresponds to a 
single plant sample. (C) Anti-PPV CP (CP) western blot assay of upper 
non-inoculated W. clevelandii leaves (22 dpi); each lane corresponds to a 
single plant sample. Ponceau red-stained blots are shown as loading 
control. 

doi:10.1371/journal.ppat.1003985.g004 

with FG-6,7-AA (FG) substitution in the IxFG motif or W63A 
E67A (WE) substitution in the FxxLE motif (Figure 3B), in 
accordance with the finding that the N-terminal region is 
dispensable for potyvirus PI processing [18]. We tested whether 
the Thosea asigna virus 2A (T2A) "self-cleaving" peptide [36] could 
overcome PI protease defects and restore Pl-HCPro separation. 
We confrrmed that T2A insertion between PI S259A and HCPro 
(ST2A) restored precursor processing, which, considering the 
barely detectable level of the predicted 48.9 kDa uncleaved 
product, was more efficient than with the wild- type PI construct. 

Involvement of PI conserved motifs in viral infectivity 

To test their effect on viral infectivity, the PI alanine 
substitutions used in in vitro translation were introduced into 
binary vectors bearing the full-length PPV cDNA. JV. clevelandii 
plants were agro-infiltrated with Agrobacterium strains harboring 
pSN-PPV, pSN-PPV API, or pSN-PPV plasmids with point 
mutations in PI, and GFP ffuorescence was monitored to follow 
the infection process of the derived viral clones. Infectious clones 
with full deletion (API) or point mutations in the PI sequence (FG 
and WE) had characteristic ring-shaped GFP foci with no 
detectable fluorescence in the center compared to the wild-type 
virus (PPV) (Figure 4A). A similar phenotype was observed in a 
Potato virus A (PVA) clone with the GFP gene inserted into the P 1 
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region [37], suggesting that in this recombinant potyvirus, PI 
function was also affected. Plants challenged with clones harboring 
the S or the HY mutations, which had no proteolytic activity in the 
in vitro translation assays, did not show GFP fluorescence, infection 
symptoms or CP accumulation (Figure 4A-C and not shown). 

Although PI VE-189,190-AA replacement led to no detectable 
cleavage in WGE experiments (Figure 3), a faint CP accumulation 
signal (compared to wild-type PPV) was visible in leaves agro- 
infiltrated with the VE clone harboring the same mutations, and 
the CP signal was gready enhanced in systemically infected leaves 
(Figure 4B and C). These leaves showed the characteristic ring- 
shaped GFP foci typical of the other infectious PI mutants 
(Figure 4A). The results led us to hypothesize that partial reversion 
might occur, restoring PI protease activity. Western blot 
performed using anti-PPV HCPro antibody confirmed correct 
Pl-HCPro processing in the upper non-inoculated leaves of plants 
agro-inoculated with the VE mutant clone (Figure S2A). Samples 
of systemi[:ally infected leaves from these plants were further 
subjected to RT-PCR amplification of a PPV genome fragment 
spanning the mutations introduced. Though the El 90 A mutation 
was maintained, we detected reversion of V189A to the wild-type 
residue at 22 dpi (Figure S2B). It is likely that, in contrast to clones 
S and HY, clone VE maintainc-d minimal PI self-processing 
activity sulficient to initiate viral replication and to select viral 
mutant progeny with improved cleavage efficiency. This is 
consistent with a report of a TEV clone in which disruption of 
the (at that time undescribed) VELI motif appeared to preclude PI 
self-processing in vitro. A virus harboring the same mutation was 
infectious, however, and able to move systemically, with a marked 
delay compared to the parental virus [20]. 

At the same time point at which PI VE- 1 89, 1 90-AA reversion 
was detected (22 dpi), the rest of the mutations introduced were 
stably maintained in both FG and WE infectious clones (verified 
by RT-PCR and sequencing, not shown), confirming that 
mutations of N-terminal motifs are less detrimental than protease 
defects. 

PI protease defects impair the silencing suppressor 
activity of HCPro 

Viabflity of the TEV clones altered by mutations in the PI 
protease can be restored by insertion of a surrogate cleavage site 
recognized by the TEV NIa protease [20], or in transgenic plants 
expressing the Pl-HCPro cistron [32]. To test whether these 
rescue strategies complement defects in HCPro function rather 
than in PI, we pc'rform(;d a transic'nt RNA silencing assay in J\f. 
benthamiana. Leaves were c:o-infiltrated with an Agrobacterium strain 
bearing p35S:GFP as a silencing reporter and Agrobacterium strains 
containing the PPV silencing suppressor HCPro preceded by the 
wild- type PI (pWT), by PI with an alanine replacement of the 
catalytic S259 (pS) or by PI S259A plus T2A (pST2A). GFP 
fluorescence was visible in all agroinfiltrated patches at 3 days 
post-agro-infiltration (dpa) (not shown). At 6 dpa, there were no 
differences between the empty control and the S construct. In 
contrast, bright fluorescence as a result of silencing suppression 
was maintained when HCPro was effectively released by PI 
protease activity or by the ribosome skipping mechanism of T2A 
(Figure 5 A). 

PPV viral clone ST2A, into which we inserted the T2A peptide 
sequence between PI S259A and HCPro, was able to initiate viral 
replication (6 dpi; Figure 5B) and move systemically despite the PI 
protease-inactivating mutation (Figure S2C; maintenance of PI 
S259A substitution was verified by RT-PCR and sequencing, not 
shown). The presence of T2A was nevertheless insufiicient to fuUy 
complement viral defects, since the ST2A GFP fluorescence 



phenotype and CP accumulation levels differed considerably from 
wild-type PPV (Figure 5B, Figure S2C and D). 

To further validate these results, A. ihaliana Col-0 and the RNA 
silencing-defective dcl2/3/4 line were agro-inoculated with the 
PPV S259A mutant (S), using wild-type PPV as control. As 
predicted, wild-type PPV infected 16 of 16 ( halk-ngcd plants of 
both the mutant line and its wild-type background, at 23 dpi. 
Although we could not detect replication of the S mutant in Col-0 
by GFP fluorescence or western blot analysis, GFP fluorescence 
was observed in 11 of 16 dcl2/3/4 plants (Figure 5C). In spite of 
the cleavage-disturbing mutation in the S cDNA clone, Pl-HCPro 
proteolytic separation was rescued in dcl2/3/4 plants inoculated 
with PPV S259A, as verified by western blot assay of HCPro 
(Figure 5D). RT-PCR analysis of samples collected at 23 dpi 
confirmed that the original serine and the protease activity were 
restored in the viral progeny, since the wild-type serine 259 codon 
AGC, which was mutated to alanine GCC in the PPV S259A 
cDNA clone, was further mutated to serine UCC (Figure 5E). The 
fact that viral reversion mutations are promptiy selected in PPV 
S259A-infected plants suggests that PI, in addition to the release of 
an active silencing suppressor, has further function(s). 

The first 164 N-terminal residues are dispensable for PPV 
PI self-processing and act as a negative regulator in in 
vitro cleavage assay 

According to MEROPS classifK:ati()n [38], potyviral PI serine 
protease belongs to subclan PA(S), whose representative member, 
trypsin, is synthesized as the inactive precursor trypsinogen, with 
a disordered loop partially obstructing the substrate-binding cleft 
[39]. Supported by the intrinsic disorder prediction of the PI N- 
terminal rc^gion and since a non-viral factor is needed for PI 
protease' acti\'ati()n [18], we tested whether PI also fits the trypsin 
model. Pre\'ious studies on TEV [40], sequence alignment, 
secondary structure predictions, intrinsic disorder confidence, as 
weU as the fmding that in PPV, the VE- 189,1 90-AA substitution 
disturbed PI self-processing, were considered in choosing PPV PI 
residues T162 and SI 70 as truncation points for a preliminary 
trial (Figure 6A). N-terminal deletion constructs were made by 
removing the P 1 sequence upstream of the codon for each amino 
acid selected, except for the initial methionine. To test PI 
protease activity, we used in vitro translation in WGE and rabbit 
reticulocyte lysate (RRL) systems. The full-length PI (WT) 
construct released the mature 35.3 kDa PI protein in WGE but 
not in the RRL system (Figure 6B, single asterisk), as anticipated 
[18]. Nonetheless, the T162 deletion construct successfully self- 
cleaved in both WGE and RRL systems, as shown by release of 
the 16.7 kDa PI processed fragment (Figure 6B, double asterisks). 
The SI 70 construct lacked activity in both in vitro translation 
systems, and only its uncleaved 27.0 kDa precursor was 
detectable (Figure 6B), suggesting that the PI residue delimiting 
the N-terminal minimal protease domain is located between 
positions 162 and 169. 

To fine-map this boundary, additional single amino acid 
deletions were engineered, transcribed, and tested by in vitro 
translation. As a further control, the catalytic S259A mutation was 
included in the T162 construct to rule out misleading non-specific 
protein degradation. Efficient cleavage activity was maintained in 
both WGE and RRL after deletion of P 1 N-terminal amino acids 
2-164 (Figure 6C, R165 construct), confirmed by release of the 
= 16 kDa PI mature fragment. Further truncations of the protease 
domain led to the drastic disappearance of protease activity in 
RRL, and gradually decreasing efficiency in WGE. The control 
protease catalytic mutant T162-S showed only a band corre- 
sponding to the unprocessed 27.8 kDa precursor. The results 



PLOS Pattiogens | www.plospatfiogens.org 



5 



Marcfi 2014 | Volume 10 | fssue 3 | el 003985 



PI Protease Modulates Host Defense Responses 



A 

GFP 
pS 

pST2A| 
pWT 



GFP 



Pl*|HCProh - 
Pl*|HCProh 




pST2A 



B 



P1 HCProh 



> 

X CL 
D. 



Col-0 dcl2/3/4 




^ > 

H Q- 
OT CL 



CP 




259 

wtPPV ATGAGCGGA 
L S J 




ATGgcCGGA 

L A J 




ATG_tcCGGA 

L S J 



Figure 5. RNA silencing suppression is impaired by the iacl< of 
separation between PI and HCPro. (A) The transient RNA silencing 
assay was done by co-infiltrating W. benthamiana leaf tissue with an 
Agrobacterium p35S:GFP culture and cultures of Agrobacterium without 
plasmid (fD) or containing pSN.5 P1-S (pS, producing PI S259A+HCPro), 
pSN.5 P1-ST2A (pST2A, producing PI S259A+T2A+HCPro), or pSN.5 PPV 
(pWT, producing wild-type Pl4-HCPro). Picture was taken under UV light 
(6 dpa). Scale bar, 1 cm. (B) N. develandii plants were agro-inoculated 
with wild-type PPV and viral clones S (PI S259 replaced by alanine) and 
ST2A (PI S259A-FT2A peptide). Viral accumulation in inoculated leaves 
was assessed by anti-PPV CP (CP) immunoblot analysis (6 dpi); each 
lane represents a pool of samples used for signal quantification. Relative 
CP signal intensities were plotted using average PPV value equal to 1 00. 
Histograms show mean ± SD (n = 4 samples/condition, from two 
independent Agrobacterium cultures); ** p<0.01, Student's f-test. 
Ponceau red-stained blot is shown as loading control. (C) Symptoms 
and GFP fluorescence of Arabidopsis dd2/3/4 plant agro-inoculated with 
PPV mutant clone S (PI S259 replaced by alanine). Pictures were taken 
in an epifluorescence microscope at 23 dpi. Scale bar, 1 cm. (D) Anti- 
PPV HCPro (HC) western blot assay of agro-inoculated Arabidopsis 
plants. Samples of upper non-inoculated leaves were collected (23 dpi). 
H, healthy dcl2/3/4 plant; PPV, pool of 16 Col-0 plants challenged with 
wild-type PPV; S, plants challenged with viral clone S, where Col-0 
samples are pools of four plants and dd2/3/4 are pools of two/three 
dd2/3/4 plants showing GFP fluorescence. Ponceau red-stained blot is 
shown as loading control. (E) RT-PCR of A. thaliana dd2/3/4 plants 
challenged with viral clone S (23 dpi). A PPV cDNA fragment 
encompassing the mutated regions (nt 1-1197 of the PPV genome) 
was amplified by RT-PCR from a pool of three plants and analyzed by 
DNA sequencing. Sanger sequencing results of the RT-PCR fragment 
and of the inoculated plasmid DNA are shown. The wild-type PPV 
sequence is shown as reference. Nucleotides that differ from the wild- 
type sequence are indicated in lower case. The nucleotide that leads to 
reversion of the amino acid codon is underlined. Encoded amino acids 
are indicated beneath each nucleotide sequence and between box- 
drawing characters. The mutated residue is in italic (alanine), reverted 
residue underlined. 

doi:1 0.1 371 /journal.ppat.1 003985.g005 
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Figure 6. Mapping of PI minimal protease domain and 
assessment of cleavage activity in eukaryotic translation 
systems. (A) Diagram of the DNA constructs used in the transcription 
reactions, as detailed in Figure 3, and partial alignment of PI from PPV 
and three reference potyviruses (see GenBank accession numbers in 
Table SI, in Text SI). Beside each virus name, PI size is indicated. 
Positions of the first and last PI amino acid shown in the alignment are 
indicated. Numbers above the alignment indicate the position of PPV 
PI residues selected for a preliminary trial. The star marks VI 64, the first 
amino acid maintained in the mutant clone used for the viral replication 
assay. The VELI motif is boxed. Residues aligning to the PPV PI minimal 
protease domain (see below) are shown in black letters, the remaining 
residues in gray letters. Alignment conservation bars are displayed. 
JNet, PsiPred and SSPro secondary structure predictions are presented, 
a-helices, black ovals; p-sheets, grey arrows. PPV DP2, disorder 
prediction of PPV PI sequence according to DIS0PRED2, where "+" 
are disordered and "-" ordered residues. PPV DP2_conf DIS0PRED2 
prediction confidence. (B) PI self-cleavage was evaluated by in vitro 
translation. N-terminal deletions were named according the first amino 
acid maintained (in addition to the initial methionine) in the truncated 
PI protein tested. Labeled translation products were resolved by 
tricine-SDS-PAGE and the ^^S signal detected. Products of reactions 
performed using wheat germ extract (WGE; top panel) and rabbit 
reticulocyte lysate (RRL; bottom). Molecular weight markers are shown 
(left). Double asterisks indicate the 16.7 kDa PI fragment released by 
the T162 construct. Single asterisk indicates the full-length PI protein 
(35.3 kDa) released by the WT construct. (C) Fine mapping of the PPV 
PI protease domain. Nomenclature as in (B). T162-S also presents 
alanine replacement of the catalytic S259 residue. WGE reaction, top 
panel; RRL, bottom. Star marks VI 64, as the corresponding PI N- 
terminal deletion was introduced into a PPV cDNA clone for the viral 
replication assay. Molecular weight markers are indicated (left). 
doi:10.1371/journal.ppat.1003985.g006 



indicate that the protease catalytic domain is correctly folded even 
in RRL, and that the first 164 N-terminal residues are not only 
dispensable for this activity, but show an antagonistic effect on PI 
self-processing in RRL. 
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Removal of the PI N-terminus enhances early 
amplification of PPV 

The PPV cDNA clone pSN-PPV PI Pro [VI 64], lacking PI 
amino acids 2-163 (PlPro, Figure 7A), was engineered to evaluate 
how a PI protease free of its antagonistic N-terminal region affects 
viral replication. We deleted PI residues upstream of position 164, 
since the in vitro cleavage assays showed that the VI 64 construct 
was the most efficient PI truncation tested (Figure 6C, star). 
Various reporter genes cloned into viral infectious cDNAs were 
used to quantify potyviral genome replication rates [41,42]. 
Taking advantage of the GFP marker inserted in our PPV clones, 
fluorescence intensity (FI) signals were used to monitor viral 
amplification kinetics. Compared to the wild-type PPV, FI levels 
were significantly higher in leaves agro-inoculated with PlPro 
clone at 46 h post-agro-inoculation (hpi) and 54 hpi (3.6- and 3.1- 
fold, respectively). At later time points, while the FI of wild-type 
PPV continued to increase over the 6-day time course, that of 
PlPro slowed at 3 dpi and dropped to 0.4 times the level of 
parental PPV at 148 hpi (Figure 7B). To support the GFP FI 
results, growth dynamics were analyzed in viral RNA by RT- 
qPCR (Figure 7C and D). Fluorescence quantification values 
correlated positively with the amounts of PPV (+)RNA and 
(-)RNA (Spearman's R.s = 0.857 and Rs = 0.976, respectively). 
Viral amount at the end of the growth curve was further assessed 
by anti-PPV CP western blot, and confirmed that significantly less 
PlPro virus accumulated than wild-type PPV at 148 hpi 
(Figure 7E). These data demonstrate that initiation of PPV 
amplification is delayed transiently in the presence of the P 1 N- 
terminal end, which is necessary to maintain higher viral 
accumulation rates in the long term. 

To test whether the PPV PlPro viral decline depends on defects 
in RNA silencing suppression, we performed a transient agro- 
infiltration assay in jV. benthamiam and N. clemlandii plants. GFP was 
used as silencing reporter, and was stably maintained in leaf 
patches co-expressing HCPro preceded either by the wild-type PI 
(pWT) or by the PI lacking amino acids 2-163 (pPlPro), as 
evaluated by GFP FI and western blot analysis (Figure S3). This 
suggests that deletion of the P 1 N-terminal region does not impair 
HCPro silencing suppressor activity. 

Removal of the PI N-terminus enhances PPV symptom 
severity and is associated with pathogenesis-related (PR) 
protein accumulation 

To confirm and complement the early growth kinetics, we agro- 
inoculated jV. clevelandii plants with PPV PlPro and evaluated 
systemic infection. Symptoms in these plants, which were soon 
much more severe than those of plants inoculated with wild-type 
PPV, included marked stunting and necrotic lesions in the center 
of the larger chlorotic spots. As for the API and PI mutants 
(Figure 4), GFP fluorescence faded with progressing focus 
expansion in P 1 Pro-infected leaves, with no detectable signal in 
the center, and showing the characteristic ring shape before the 
appearance of necrosis (Figure 8A). RT-qPCR analysis confirmed 
that in systemically infected leaves (21 dpi), despite more severe 
symptoms, less viral RNA accumulated in P 1 Pro-infected plants 
than ill wild-type PPV samples (Figure 8B). The identity of 
infecting viruses was confirmed by RT-PCR that encompassed the 
deletion, and by sequencing (not shown). 

At the protein level, PI Pro-infected plant extracts were character- 
ized by a specific band migrating at —35-37 kDa, which was 
unappreciable in wild-type PPV protein extracts and was thus 
analyzed by mass spectrometry (MS) (Figure 8C). In the MALDI- 
TOF retrieved spectrum, the P 1 Pro band showed six prominent peaks 
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Figure 7. Early amplification dynamics for wild-type PPV and 
PlPro viral clones. (A) Leaves of ttiree-weel<-old N. clevelandii plants 
were agro-inoculated witfi pSN-PPV (wild-type PPV) or pSN-PPV 
P1Pro[V164] (PlPro, in wfiich PI amino acids 2-163 were deleted). 
Tfie reporter sGFP(S65T) [103] gene is present between Nib and CP 
coding sequences of both viral clones. Boxes with diagonal lines 
indicate P3N-PIP0 protein. (B) GFP fluorescent intensity (FI) from 
inoculated leaves was quantified in a fluorometer and plotted using 
average wild-type PPV value at 148 hpi equal to 100. Line graph shows 
mean ± SD (n = 16 samples/condition, from two independent Agro- 
bacterium cultures); colors as in (A). (C) Amount of viral (+)RNA from 
inoculated leaves was quantified by RT-qPCR and plotted using average 
wild-type PPV value at 148 hpi equal to 100. Line graph shows mean ± 
SD (n = 3/4 biological replicates, from two independent Agrobacterium 
cultures); colors as in (A). (D) Amount of viral (-)RNA from inoculated 
leaves was quantified by RT-qPCR and plotted using average wild-type 
PPV value at 148 hpi equal to 100. Line graph shows mean ± SD (n = 4 
biological replicates, from two independent Agrobacterium cultures); 
colors as in (A). (E) Anti-PPV CP (CP) western blot assay of inoculated 
leaves (148 hpi); each lane represents a pool of samples used for signal 
quantification. S, sample of plants inoculated with the control viral 
clone pSN-PPV P1-S. Relative CP signal intensities were plotted using 
average PPV value equal to 100. Histograms show mean ± SD 
(n = 4 samples/condition, from two independent Agrobacterium cul- 
tures); *** p<0.001. Student's f-test. Ponceau red-stained blot is shown 
as loading control. 

doi:10.1371/journal.ppat.1003985.g007 



(Figure 8C), which were analyzed by MS/MS to obtain their peptide 
fragmentation fmgerprint. Database searching allowed assignment of 
the m/z 1203.59 fragment to residues 326-335 of M tahacum acidic 
PR-2 isoform GI9 (GenBank accession no. P23547.1). After removal 
of the signal peptide, this protein has a reported molecular weight of 
34.8 kDa [43], in good agreement with the electrophoretic mobility of 
the PI Pro-specific band. Sequence of the other five major peaks was 
defined de novo (Figure S4A), since database identification of was 
unsuccessfiil due to unavailability of the host genome sequence 
(jV! clevelandii). Based on MS/MS spectra, we verified that all five 
matched a glucan endo-l,3-(5-D-glucosidase, EC 3.2.1.39, belonging 
to class II of the PR-2 family [44,45] , and that the minor amino acid 
changes identified were consistent with acidic PR-2 sequence 
variability in Solanaceae species (Figure S4B). 
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Figure 8. Symptoms and viral accumulation of PPV lacking P1 
N-terminal residues. N. clevelandii plants were agro-inoculated with 
wild-type PPV or PI Pro, pSN-PPV P1Pro[V164] in which PI amino acids 
2-163 were deleted; data were collected at 21 dpi. (A) Pictures of N. 
clevelandii systemic leaves. Boxed leaf areas were magnified and their 
GFP fluorescence, acquired by laser scanning, is shown. Scale bar, 
2.5 cm in whole leaf pictures, 0.4 cm in details. (B) Amount of viral 
(+)RNA from systemically infected leaves was quantified by RT-qPCR 
and plotted using average wild-type PPV value equal to 100. 
Histograms show mean ± SD (n = 4 biological replicates, from two 
independent Agrobacterium cultures); *** p<0.001. Student's f-test. (C) 
Total protein extracts from systemic leaf samples were resolved by 
glycine-SDS-PAGE and Coomassie blue stained, molecular weight 
markers are indicated (left). The =35-37 kDa PI Pro-specific band was 
analyzed by MALDI-TOF and its peptide mass fingerprint is shown. 
Intensity in the y axis; m/z values are indicated for the six PI Pro sample 
peptides further analyzed by MS/MS. (D) Anti-PPV CP (CP), anti-PR-2 
(PR-2) and anti-PR-3 (PR-3) western blot assays of upper non-inoculated 
leaves; each lane represents a pool of samples used for signal 
quantification. S, sample of plants inoculated with the control viral 
clone pSN-PPV P1-S. Ponceau red-stained blot as loading control of 
protein extracts. (E) Relative PR-2 and PR-3 signal intensities were 
plotted using average PI Pro value equal to 1 00. Histograms show mean 
± SD (n = 4 samples/condition, from two independent Agrobacterium 
cultures); *** p<0.001. Student's f-test. 
doi:10.1371/journal.ppat.1003985.g008 

PR proteins of different families overaccumulate in tobacco 
plants tiiat show hypersensitivity to Tobacco mosaic virus (TMV) [46- 
49] , but also after bacterial and fungal infection and in response to 
abiotic stress [50] . We first confirmed the MS results by western blot 
analysis for PR-2; next, to test whether PR-2 induction is a specific 
PI -associated defense mechanism or part of a broader stress 
response, we assessed class II PR-3 protein expression in infected jV. 
clevelandii plants (Figure 8D). We found that, despite the lower viral 
CP levels, PR-2 and PR-3 accumulation was significantly higher in 
PlPro- compared to wild-type PPV-infected plants (Figure 8E). 



Immune response to PPV PI Pro is attenuated by 
downregulation of salicylic acid signaling 

In tobacco plants, exogenous salicylic acid (SA) treatment 
induces class II PR-2 and class II PR-3 transcription [48,51,52], 
and activates the resistance responses associated to TMV-induced 
hypersensitivity [53,54]. Downregulation of SA accumulation and 
systemic acquired resistance was reported in tran.sgenic plants that 
express the bacterial salicylate-hydroxylase gene nahG [55]. 

In N. benthamiana, while wild-type PPV-infected plants appeared 
almost symptomless, PPV P 1 Pro-infected plants showed extended 
chlorosis and necrotic lesions, similar to those observed in 
jV. clevelandii (Figure 9A). We therefore used transgenic 
jV. benthamiana NahG plants [56] to evaluate the SA contribution 
in the PPV PlPro host immune response. In accordance with 
previous studies [56], downregulation of SA signaling had no 
appreciable effect on the wild-type PPV phenotype. In contrast, in 
PI Pro-infected NahG plants, systemic symptom severity was 
attenuated (10 dpi; Figure 9A and B). This result was supported by 
a sharp reduction in PR-2 protein accumulation in P 1 Pro-infected 
NahG plants and a weak increase in viral load, estimated by anti- 
PPV CP western blot analysis (12 dpi; Figure 9C and D). These 
data suggest that although SA-mediated antiviral pathways have 
only a minor role in wild-type PPV infection, they take part in 
PlPro immune responses. 

Several stress response proteins are upregulated in 
PI Pro-infected NahG plants 

In NahG-expressing plants (i) PI Pro-induced chlorosis was 
more accentuated than in wild-type PPV-infected plants 
(Figure 9A), (ii) although lower than the wild-type host, PR-2 
abundance in P 1 Pro samples was significantly higher than in wild- 
type PPV, and (iii) PR-2 reduction was insufficient to fuUy restore 
PlPro viral accumulation to parental PPV levels (Figure 9D), 
These findings prompted us to further investigate PI Pro-related 
defense responses after downregulation of SA signaling. To 
identify proteins whose abundance was significandy changed in 
P 1 Pro-infected NahG plants relative to wild-type PPV-infected 
NahG plants, we performed quantitative proteomic analysis using 
isobaric tag labeling (iTRAQ) and liquid chromatography (LC)- 
MS/MS [57]. A draft sequence of the M. benthamiana genome was 
released [58,59], and a search against its predicted protein 
database enabled us to identify more than a thousand non- 
redundant proteins. Of these, 23 were considered to accumulate 
differentially in PlPro versus wild-type PPV-infected plant 
samples, since they were found in both P 1 Pro biological replicates 
#A and #B, with a false discovery rate <5% as statistical cut-off 
(Figwe lOA). Gene ontology term enrichment analysis showed 
that, according to plant symptoms, the 23 dysregulated proteins 
associated significantly with the GO term "response to stress" (GO 
ID: 6950, /)<0.0001; Figure lOB). In Figure IOC, we present a 
heat map of these quantified PlPro proteins with their average 
iTRAQ^ ratios (expressed in LOG2) relative to wild-type PPV 
biological replica #A. As predicted, PlPro biological replicates 
jj^K and t^B are grouped in the same hierarchical cluster, which 
differs from PPV biological replica t^B. In PlPro samples, several 
proteins from different PR families were more abundant than in 
wild-type PPV. These include a class II PR-2 (in accordance with 
the western blot result; Figure 9C) and other acidic members, 
which are regulated in tobacco mainly by SA [50,52], and basic 
counterparts such as the basic PR-1, whose transcription is 
effectively activated by ethylene [60-64]. Abscisic acid and 
osmotic stress are reported to induce expression of basic PR-5 
[50,65], as well as of dehydrin-like proteins [66]. The contribution 
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Figure 9. Viral accumulation and immune response induction in NahG-expressing N. benthamiana. (A) Wild-type (N.b. wt) and NahG- 
expresslng (N.b. NahG) N. benthamiana plants were agro-inoculated with wild-type PPV and PI Pro viral clones. Pictures of systemic leaves were taken 
at 10 dpi. GFP fluorescence, acquired by laser scanning, is shown. Scale bar, 4 cm. (B) Details of (A) showing symptoms of wild-type and NahG- 
expressing plants infected with the PI Pro viral clone. (C) Anti-PPV CP (CP) and anti-PR-2 (PR-2) western blot analysis of upper non-inoculated leaves 
of wild-type and NahG-expressing plants (12 dpi); each lane represents a pool of samples used for signal quantification. S, sample of plants 
inoculated with the control viral clone pSN-PPV P1-S. Below, blot stained with Ponceau red as loading control of protein extracts. (D) Relative PR-2 
and CP signal intensities were plotted using average PI Pro or wild-type PPV value equal to 100, respectively. Histograms show mean ± SD 
(n = 4 samples/condition, from two independent Agrobacterium cultures). Bars with different letters are statistically significant, p<0.05, one-way 
Anova and Tukey's HSD test. 
doi:10.1371/journal.ppat.1003985.g009 



of oxidative .stress in the SA-dependent response to PI Pro is 
underlined by the detection of a peroxidase and a catalase; the 
tobacco homologue of the latter was initially identified as a SA- 
binding protein [67]. Only three of the 23 difiFerentiaUy 
accumulated host proteins were downregulated in PI Pro; these 
include two plastocyanins and a CP12-like accession, proteins 
essential for photosynthesis [68,69] . These data suggest that in the 
transgenic host, an important component of SA signaling is 
maintained and/or that its downregulation is partially compen- 
sated by alternative defense components. 

Discussion 

Besides the self-cleavage activity intrinsic to its C-terminal end, 
other activities of potyviral PI and its real contributions to viral 
replication remained vague, so much so that PI earned the 
appellative "mysterious protein" in a recent review [24]. 

In concomitance with the discovery that HCPro inhibits 
posttranscriptional gene silencing, it was suggested that PI could 
act in conjunction with and strengthen the silencing suppressor 
activity of HCPro [25,26]. A characteristic ring-shaped GFP focus 
phenotype was related in tobamovrrus infections with deficiency in 
RNA silencing suppression [70,7 1] . The same GFP phenotype was 
distinctive of all our PPV infectious clones with PI sequence 
deletions and point mutations. In contrast, our Arabidopsis dd2/3/ 4 



complementation data demonstrate a major function for PI 
independent of RNA silencing suppression. The ring GFP 
phenotype of PI mutants might thus be related not only to a 
silencing suppression defect, but also to other altered functions. 
Recently, the silencing suppressor-enhancing effect of PI was 
attributed to cu elements that improve translation efficiency, rather 
than to complementary activity of the PI protein [72]. 

Our results nonetheless show that PI is unquestionably 
involved in RNA silencing by limiting HCPro function. 
Deficiencies in polyprotein processing at the Pl-HCPro junction 
abrogated PPV infectivity in wild-type plants, as also reported for 
TEV [20], and were partially complemented by insertion of the 
T2A "self-cleaving" peptide. To further define the nature of this 
viral weakness, we used a transient RNA silencing assay to show 
that for correct sUencing suppression, PI must be effectively 
separated from HCPro, whether mediated by its own protease 
domain or by the ribosome skipping mechanism of T2A. In 
plants with impaired antiviral RNA silencing machineries the 
PPV S259A clone, with a cleavage-inactive PI protease, 
recovered infectivity. The selection of viral revertant progeny 
demonstrated that strong evolutionary pressure for proteolytic 
competence of PI was maintained, despite the host mutant 
background. The lack of separation between PI and HCPro 
probably impairs not only the RNA silencing suppressor activity, 
but also other important viral functions. 
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Figure 10. Quantitative proteomic analysis of PPV P1 Pro- 
infected NahG plants. (A) Venn diagram of proteins whose 
abundance was significantly changed In PI Pro biological replicates 
#A and #B, relative to wild-type PPV sample #A. A false discovery rate 
of <5% was applied as statistical cut-off Numbers In the non- 
overlapping areas represent unique significant proteins of each 
biological replica; the number In the overlapping area, 23, shows the 
common significant proteins for further analysis. (B) Map of gene 
ontology terms significantly enriched In 23 dysregulated proteins 
common to both PI Pro replicates. Relevant biological process nodes 
are colored from yellow to red, where red Is the most significantly 
overrepresented category. In white, non-slgnlflcant node. (C) Heat map 
of 23 significantly differentially regulated proteins common to PI Pro 
replicates #A and #B. Fold change was calculated using PPV #A values 
as 1, and expressed In L0G2. Quantification values of PPV biological 
replica #B are shown as control. Sample tree was built by hierarchical 
clustering analysis. Beside each W. benthamiana accession ID, protein 
Isoelectric point Is Indicated. Brief protein descriptions are reported on 
the right; PR proteins are grouped by families according van Loon ef al. 
[50]. 

dol:10.1371/journal.ppat.l003985.g010 



Of the mature proteins encoded by the potyviral genome, PI 
presents the greatest variability in length and in amino acid 
sequence [73]. In silico analysis showed that the limited primary 
sequence conservation of the P 1 N-terminal region is associated 
with residues predicted to be part of intrinsically disordered loops. 
The relevance of unfolded regions in the proteolytic maturation of 
viral polyproteins has been shown [6,7], and flexible loops in 
peptidase precursors act in many cases as protease activation 



switches [39]. Previous results and those reported here show 
successful self-cleavage activity of potyviral PI after in vitro 
translation only in the WGE system, but not in RRL. In 
consequence, it was suggested that a non-viral factor, present in 
WGE and absent in RRL, might be involved in correct folding of 
the protease domain or presentation of the cleavage site [18]. Our 
evaluation of the cleavage performance of PI N-terminal end 
truncations demonstrates that the PI minimal protease domain is 
active in WGE but also in RRL, suggesting that both eukaryotic 
translation systems support correct folding and activity of the P 1 
protease domain. Hence, the PI N-terminal residues are not only 
dispensable for HCPro release, but also act as a negative regulator 
of P 1 self-processing. 

In the context of viral infection, removal of the PI protease 
antagonistic extension accelerated early viral replication and was 
followed by enhancement of symptom severity, with the appear- 
ance of plant stunting and necrotic lesions that did not 
characterize wild-type PPV. The greater aggressiveness did not 
parallel viral loads, however, as viral RNA and CP levels in plants 
systemically infected with the deletion mutant were lower than in 
wild-type PPV-infected plants. Lack of positive correlation 
between symptom severity and fitness is reported for several viral 
systems, including potyviral infections [74—76]. This can be 
justified by the crossing of a virulence threshold, which would 
result in higher induction of immune responses or excessive host 
debilitation [77,78]. Accordingly, the necrotic phenotype that 
characterizes PPV P 1 Pro was associated with overaccumulation of 
class II PR-2 and PR-3 proteins. Acidic PR-2 and PR-3 family 
members are indicators of SA-mediated responses to biotic stress 
and are hallmarks of systemic acquired resistance [48,50,79]. In a 
transgenic NahG-expressing plants, a decrease in SA signaling 
attenuated PPV PI Pro symptom severity. PI Pro-induced chlorosis 
was nonetheless more pronounced than in wild-type PPV-infected 
plants. Complex hormone crosstalk follows plant infection [80], 
and compounds other than SA operate in the defense against 
pathogens [81-83]. Quantitative proteomic analysis of NahG 
plants infected with the PI Pro viral clone allowed us to identify a 
large set of stress response-associated proteins whose abundance 
was significandy altered compared to wild-type PPV. This includes 
upregulation of members of different PR families whose 
transcription is reported to be activated by SA but also by 
ethylene, and further modulated by abscisic acid [50], as well as 
rearrangement of the antioxidant system, and downregulation of 
photosynthetic components, similar to other studies [84-86]. 
Coordinated action of both SA-dependent and -independent 
pathways probably contributes to the immune response activated 
by PPV PI Pro. In turn, the N-terminal extremity of the PI 
protease, absent in PI Pro but maintained in wild- type PPV, helps 
to bypass induction of host defense responses. Several plant viruses 
interfere with salicylate pathways [87-89], and it would be of 
interest to determine whether PI also has an active role in 
suppressing basal host defenses independent of its HCPro activity- 
modulating effect. 

Considering that (i) defects in PI self-cleavage preclude viral 
viability, (ii) viral RNA silencing suppression is impaired by the 
lack of separation between PI and HCPro, (iii) PI residues 1-164 
are predicted to be mainly disordered, and negatively affect Pl- 
HCPro processing in RRL, (iv) PPV early amplification dynamics 
are enhanced by removal of the P 1 N-terminal end, (v) increased 
defense responses are associated with deletion of the PI N- 
terminus, and (vi) removal of the PI N-terminus reduces viral 
accumulation, we propose a model by which PI acts to fine-tune 
potyviral replication by sensing specific host effectors. The 
presence of these cofactors leads to activation of the PI protease 
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domain and release of the silencing suppressor HCPro. The virus 
can thus efTectively counteract host antiviral RNA silencing 
defenses and replicate successfully, but it can also stimulate 
additional defense responses. In an alternative scenario in which 
the cofactor is limited, PI -HCPro separation is restricted, RNA 
silencing suppression activity is disturbed, and some decrease in 
viral replication can therefore be predicted. We speculate that 
this mechanism holds viral replication below a level that would be 
detrimental to the host cell and reduces triggering of host 
immune responses, thus maintaining higher long-term replicative 
capacity. 

An example of restriction in a virus infection step to avoid host 
defense responses is the regulation of autoprocessing of the NS2 
protein of the pestivirus Bovine viral diarrhea virus (BVDV), which is 
dependent on the host chapcrone Jiv [90] . Limiting amounts of Jiv 
cofactor result in accumulation of uncleaved NS2-3 and RNA 
replication arrest. Like the necrogenic interaction of PPV P 1 Pro, 
the disengagement of NS2 self-cleavage from cellular Jiv leads to a 
switch from a non-cytopathogenic to a cytopathogenic biotype, 
underlining the importance of the temporal modulation of NS2-3 
processing [91,92]. 

To understand the regulation of PI -HCPro processing, 
identification of the PI activating host effector will be a priority 
for future studies. Moreover, PI protease multifunctionality was 
postulated [20,93] and indications that leader proteases retain 
activities independent of proteolytic processing can be found in 
several viruses [94-96] . Given the limited infection efficiency of 
the PPV API viral mutant and the effects of mutagenesis of PI 
conserved motifs, it is likely that PI protein or its RNA 
coding sequence has additional unrevealed roles in potyviral 
replication. 

Materials and Methods 

Sequence analysis 

Potyviral PI sequences were visualized and analyzed using 
Jalview [97]. Multiple sequence alignments were made using 
MUSCLE [98]. Protein disorder predictions were made with 
DISOPRED2 [33], setting a false positive rate threshold of 10%, 
and MetaDisorderMD2 [34]. Secondary structure predictions 
used the JNet algorithm [99] included in the Jalview package, 
PsiPred [100] and SSpro [101]. 

Plant nnaterials 

Micotiana clemlandii, N. benthamiana, N. benthamiana NahG [56] 
(kindly provided by Prof H.S. Guo, Chinese Academy of Sciences, 
Beijing, China, and Dr. F. Tenllado, CIB-CSIC, Madrid, Spain), 
Arabidopsis thaliana Col-0 and its triple mutant line dcl2-l dcB-l 
dcl4-2 [102] (kindly provided by Prof. J.C. Carrington, CGRB, 
Oregon State University, CorvaUis, OR, USA) were used. Mcotiana 
plants were grown in a greenhouse maintained at a 16 h light/8 h 
dark photoperiod, temperature range 19-23°C. Arabidopsis plants 
were grown in the greenhouse or in vitro (see conditions below). 

DNA plasmids and constructs 

A full-length cDNA copy of a PPV isolate [30], tagged with 
sGFP(S65T) [103] and inserted in the pBINPPV-NK-GFP binary 
plasmid, was reported [104]. Cloning in binary vectors is often 
constrained by the availability of unique restriction sites or, in the 
case of the Gateway technology [105], the maintenance of X 
integrase recombination sites in the final constructs. To overcome 
these limitations, pBINPPV-NK-GFP was engineered to obtain 
the pSN-ccdB plasmid. Once linearized, pSN-ccdB is suitable for 
one-step seamless replacement of the PPV 5'UTR and the PI 



cDNA sequence by Gibson assembly [106]; it was therefore used 
as a backbone for aU the viral cDNA clones in this study. A 
detailed description of the viral vectors, together with in vitro 
translation and agro-infiltration constructs, can be found in Text 
SI of the Supporting Information. 

In vitro transcription and translation 

In vitro transcription reactions were performed with the TV- 
Scribe Standard RNA IVT kit (CELLSCRIPT), including final 
DNase I digestion. RNA was purified by organic extraction/ 
ammonium acetate precipitation. Quality was assessed by 
NanoDrop (Thermo Fisher Scientific) and by gel electrophoresis, 
with final concentration adjusted to 1 |j,g/|J,L. In vitro translation 
was carried out in the presence of L-['^'S] methionine and L- 
[ 'S] cysteine (PerkinElmer) using the WGE or the RRL transla- 
tion system (Promega), according to manufacturer's instructions. 
Samples were resolved in 12% tricine-SDS-PAGE [107] and the 
signal detected by Phospholmager or autoradiography. 

Arabidopsis in vitro agro-infection 

Arabidopsis seeds were surface steriUzed, kept at 4°C for three 
days, and sown in vitro on Murashige and Skoog medium [108] 
with MES and vitamins (Duchefa), adjusted to pH 5.7 with KOH 
and supplemented with 1% (w/v) sucrose and 0.7% (w/v) Bacto 
Agar (Difco). Six days after germination, plants were transferred to 
half-strength Murashige and Skoog m(;dium with MES and 
vitamins, adjusted to pH 5.7 with KOH and supplemented with 
1% (w/v) sucrose and 0.8% (w/v) Bacto Agar. At 9-10 days post- 
germination, agro-infection was performed as follows. Agrobacterium 
cultures were induced as described [28], the final ODbqo was 
adjust("d to 1.5, forceps were dipped into the culture, and the two 
youngest plantlet leaves were pierced. Forceps were flame- 
sterilized between cultures of different bacterial clones. Plates 
were sealed with gas-permeable tape (MiUipore) and maintained in 
a growth chamber at 16 h light/8 h dark photoperiod, temper- 
ature 21±1°C, 60% relative humidity. 

Nicotlana plant agro-infiltration 

The transient RNA silencing assay and agro-infiltration of 
jV. benthamiana and N. clevelandii plants with A. tumefaciens strain 
C58C1 bearing the indicated plasmids, were as described [28]. 
The viral replication assay was conducted in three-week-old JV. 
clevelandii plants following reported agroinfiltration and sampling 
guideUnes [109], with the exception that a saturating concentra- 
tion oi Agrobacterium was used. 

GFP fluorescence imaging and quantification 

GFP fluorescence was monitored under long-wavelength UV 
hght (Black Ray, model B 100 AP) and photographed using a 
Nikon DIX digital camera with a 62E 022 filter. Alternatively, 

images were acrjuired under a MZ FLIII epifluorescence 
microscope (Leica) using a GFP3 filter (excitation and barrier 
filters at 470/40 nm and 525/50 nm, respectively) and photo- 
graphed with an Olympus DP70 digital camera. When indicated, 
GFP fluorescence was acquired by laser scanning (Typhoon 9400, 
GE Healthcare; laser 488 nm, intensity 450 V; 526 nm short-pass 
emission filter). 

GFP fluorescent intensity quantification was carried out placing 
individual 5.0 mm-cUameter leaf discs in a black 96-well plate. 
Signals were acquired in an Appliskan (Thermo Fisher Scientific) 
or Victor X2 (PerkinElmer) plate readers with the following 
settings: measurement time 500 ms, excitation and emission 
wavelengths of 485/10 nm and 535/20 nm, respectively. 
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Western blot assays 

Plant tissue was ground in a mortar in liquid nitrogen or 
homogenized in a TissueLyzer bead mill (Qiagen). Total proteins 
were extracted in 150 mM Tris-HCl pH 7.5, 6 M urea, 2% (w/v) 
SDS, 5% (v/v) glycerol and 5% (v/v) P-mercaptoethanol; heat 
denatured (96°C, 5 min) and centrifuged (14000 rpm, 10 min) to 
remove cell debris. Protein samples were separated by glycine- 
SDS-PAGE and electroblotted onto a nitrocellulose membrane. 
Ponceau red staining was used to control protein loading 
equivalence. Proteins were detected using anti-PPV CP and anti- 
PPV HCPro rabbit sera, and anti-GFP monoclonal antibody 
(clones 7.1 and 13.1, Roche) as primary antibodies. Antibodies 
raised against tobacco class II PR-2 [44] and class II PR-3 proteins 
[47] were kindly provided by Dr. T. Heitz (IBMP-CNRS, 
Strasbourg, France). Horseradish peroxidase-conjugated goat 
anti-rabbit IgG (Jackson) or sheep anti-mouse IgG (GE Health- 
care) were used as secondary- antibody. Immunostained proteins 
were visualized by enhanced chemiluminescence detection with a 
LiteABlot kit (Euroclone). For signal quantification, chemiolumi- 
nescence was acquired in a ChemiDoc XRS imager (BioRad) and 
analyzed with ImageJ software [110]. 

RT-PCR and RT-qPCR 

Total RNA was extracted with the FavorPrep Plant Total RNA 
Mini kit (Favorgen). Fragments spanning the PPV 5 'UTR, and the 
coding sequences of P 1 and the HCPro N-terminus were amplified 
with the Titan One Tube RT-PCR kit (Roche) using primers 
1595_F/1597_R (Table S2, in Text SI). When indicated, 
fragments were purified using the FavorPrep Gel/PCR Purifica- 
tion kit (Favorgen) and DNA was sequenced. 

Strand-specific quantification of PPV RNA was done for at least 
three biological rephcates per condition using tagged cDNA 
primers in the RT step [111-113], and will be detailed elsewhere 
(in preparation). Briefly, equal amounts of DNAsel-treated total 
RNA were used for cDNA synthesis using Superscript III 
(Invitrogen) and primer Q26_R or Q_29_F to transcribe cDNA 
from positive and negative PPV genomes, respectively. Technical 
triplicate qPCR reactions were prepared using HOT FIREPol 
EvaGreen qPCR Mix Plus (Soils BioDyne) in 384-well optical 
plates and run in a 790()HT Fast Real-Time PGR System (Applied 
Biosystems). Primer pairs Q27_F/Q28_R and Q30_F/Q31_R 
were used for positive and negative genome quantifications, 
respectively. The amount of target RNA in the analyzed samples 
was estimated by absolute quantification using an external DNA 
standard curve [114]. 

MALDI MS/MS and data analysis 

Plant total protein extracts were separated by glycine-SDS- 
PAGE and stained with Coomassie blue. Gel bands of interest 
were excised manually, reduced, alkylated and in-gel digested. 
The tryptic-eluted peptides were subjected to MALDI-TOF/ 
TOE analysis. Data were automatically acquired in an ABi 
4800 MALDI TOE/TOE mass specti-ometer (AB Sciex) 
and searched against NCBI non-redundant protein database, 
NCBInr_20121 1 16. Mass tolerance for precursors was set to 
±50 ppm and for MS/MS fragment ions to ±0.3 Da. The 
confidence interval for protein identification was set to ^95% 
(/)<0.05) and only peptides with an individual ion score above the 
identity threshold were considered correcdy identified. 

Manual peptide de novo sequencing was performed according 
Ma and Johnson [115]. Residues with near/isobaric masses were 
bona fide assigned according alignment consensus. Confidence of 
retrieved results was further tested with Peaks Studio software 
(BSI) [116]. 



ITRAQ sample labeling, LC-IVIS/IVIS and data analysis 

hmthamiann NahG plants were agro-inoculated with two 
independent bact('ria culture's per each viral cDNA clone (n = 2 
biological replicates); upper non-inoculated leaves were collected 
from 8 plants per culture (12 dpi). Total protein extracts were 
prepared as described for western blot assays and further purified 
by methanol/chloroform precipitation. Protein pellets were 
resuspended in 6 M guanidine hydrochloride and 100 mM 
HEPES, pH 7.5, and concentration was determined by RC DC 
assay (BioRad). Equal amounts of protein for each condition were 
trypsin-digested and labeled with iTRAQ_ Reagent Multi-plex kit 
(AB Sciex). Tags 114 and 116 were used for PI Pro biological 
replicates and 115 and 1 17 for wild-type PPV biological replicates. 
Labeled peptide samples were combined and subjected to LC- 
MS/MS analysis (three technical replicates) using a nano liquid 
chromatography system (Eksigent Technologies) coupled to a 
Triple TOE 5600 mass spectrometer (AB Sciex). MS and MS/MS 
data were processed using Analyst TF 1.5.1 software (AB Sciex) 
and searched using the Mascot Server v. 2.3.02 (Matrix Science) 
against a customized database represented by N. benthamiana 
genome-predicted proteins (available at Sol Genomics Network 
[59]) plus their corresponding reversed entries. Peptide mass 
tolerance was set to ±20 ppm for precursors and 0.05 Da for 
fragment masses. The confidence interval for protein identification 
was set to S95% (/)<0.05), and only peptides with an individual 
ion score above the identity threshold were maintained in 
quantification analysis. Proteins were considered diflferentially 
expressed if they had at least two quantified peptides and they 
were present in both PI Pro biological replicates with a false 
discovery rate <5%. LOG2 ratios relative to PPV biological 
replica #A (iTRAQ_tag 115) were visualized by MultiExperiment 
Viewer [117]. Hierarchical clustering analysis with Pearson 
correlation distance metric was used to build the sample tree 
[118]. Sequences of reliably quantified jVC benthamiana proteins 
were used as query in a WU-BLAST2 search against the TAIRIO 
protein dataset, to generate a list of homologous A. thaliana gene 
IDs. This was used as input in gene ontology term enrichment 
analysis using BinGO [119]. 

Accession numbers 

The following GenBank (http:/ /www.ncbLrJm.nih.gov) accessions 
were used in the viral sequence analysis: YMV (YP_022752.1); PW 

(NP_734369.1); PVY (NP_734243.1); WY (YP_001931974.1); 
CaYSV (YP_003208051.1); JGMV (NP_734408.1); MDMV 
(NP_734143.1); SCMV (NP_734133.1); SrMV (CAC84438. 1); 
PVA (NP_734359.1); TEV (NP_734207.1); ScaMV (NP_734123.1); 
TuMV (BAC02772.1). Plant GenBank accessions considered: N. 
tabacum class II PR-2 (P23547.1); M tahaeum class I PR-2 
(AAA63541.1); S. tuberosum class II PR-2 (CAE52322.1); S. lycopersicum 
class II PR-2 (NP_001234798.1). benthamiana accessions can be 
found at Sol Genomics Network (ftp://ftp.solgenomics.net/ 
genomes/Nicotiana_benfhamiana/annotation/). 

Supporting Information 

Figure SI Alignment and protein disorder prediction of 
the six smallest known potyviral PI sequences. GenBank 
accession numbers are reported in Table SI, in Text SI; for PPV 
PI, only resides 164—308 were considered. Amino acid back- 
ground is assigned according the ClustaLX color scheme [120]. 
Residues aligning to PPV PI minimal protease domain are shown 
in black letters and bright-colored background. The VELI motif is 
boxed; The protease conserved RG dipeptide is marked with black 
bar and the catalytic triad His, Asp and Ser is marked with stars. 
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MD2 lines show protein disorder prediction of PI sequence 
according MetadisorderMD2, where "+" are disordered and " — " 
ordered residues. DISOPRED2 prediction confidence is plotted: 
ScaMV, green line; SrMV, purple; SCMV, orange; JGMV, blue; 
MDMV, turquoise; CaYSV, magenta. 
(TIF) 

Figure S2 Pseudo-reversions of PI VE-189,190-AA mu- 
tations and infectivity of PPV ST2A viral clone. N. clevelandii 

plants were agro-inoculated with wild-type PPV and viral clone 
VE (PI VE-189,190-AA), clone S (PI S259 replaced by alanine) 
and clone ST2A (PI S259A+T2A peptide). (A) Anti-PPV HCPro 
(HC) western blot assay of upper non-inoculated leaves, collected 
at 22 dpi. For each clone, both samples shown in Figure 4C were 
analyzed together. Ponceau red-stained blot is shown as loading 
control. (B) RT-PCR analysis of plants challenged with VE 
mutant clone (22 dpi). A PPV cDNA fragment encompassing the 
mutated region (nt 1-1197 of the PPV genome) was amplified by 
RT-PCR from a pool of two plants and analyzed by DNA 
sequencing. Sanger sequencing results of the RT-PCR fragment 
and of the inoculated plasmid DNA are shown. The wUd-type 
PPV sequence is shown as reference. Nucleotides that differ from 
the wild-type sequence are indicated in lower case. The nucleotide 
that leads to reversion of the amino acid codon is underlined. 
Encoded amino acids are indicated beneath each nucleotide 
sequence and between box-drawing characters. Mutated residues 
are in italic (alanine), reverted residue underlined. (C) GFP- 
fluorescence pictures of upper non-inoculated N. clevelandii leaves 
were taken in an epifluorescence microscope (22 dpi). Scale bar, 
1 cm. (D) Anti-PPV CP (CP) western blot assay of upper non- 
inoculated leaves (22 dpi); each lane corresponds to a single plant 
sample. Ponceau red-stained blot is shown as loading control. 
(TIF) 

Figure S3 HCPro silencing suppression activity is 
maintained in the absence of the PI N-terminal region. 

Transi(-nt RNA silentdng assay was done by co-infiltrating leaf 
tissue with an Agrohaderium p35S:GFP culture and cultures of 
Agrobactenum pSN.5 Pl-S (pS, producing PI S259A+HCPro), 
pSN.5 PPV ^WT, producing wild-type Pl-hHCPro), or pSN.5 
PlPro (pPlPro, producing PI lacking amino acids 2—163-1- 
HCPro). (A) Picture of jV. benihamiana agro-infiltrated leaf was 
taken in an epifluorescence microscope (6 dpa). (B) GFP 
fluorescent intensity (FI) from J\\ henthamiana infiltrated leaf patches 
was quantified in a fluorometer and plotted using average pWT 
value at 6 dpa equal to 100. Line graph shows mean ± SD 
(«= 16 samples/condition, from two independent Agrobactenum 
cultures); colors as in (A). For each time point, the mean GFP 
FI value (n = 16) from untreated leaves was used as blank. 
Analyzed time points (dpa) in the x axis. (C) Western blot detection 
of GFP protein in JV. clevelandii samples (6 dpa); each lane 
represents a pool of samples from 2/3 infiltrated patches. N, 
non-infiltrated leaf sample. Ponceau red staining is shown as 
loading control. (D) Relative signal quantification for the 
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immunoblot shown in (C) using average pWT value equal to 
100. Histogram shows mean ± SD (n = 6 samples/condition, from 
two independent Agrobactenum cultures); ns, not significant by 
Student's Hest. 
(TIF) 

Figure S4 Nicotiana clevelandii class II PR-2 peptide 

sequences assigned by MS/MS analysis. (A) Identified 
sequence of MALDI-TOF peptides labeled in Figure 8C. Mass 
measurement accuracy of the retrieved peptides is indicated in 
parts per million (ppm). In de novo defined peptides, residues that 
differ from tobacco class II PR-2 isoform GI9 sequence (GenBank 
accession no. P23547.1) are in bold and underlined. Amino acids 
with near-/isobaric masses were bona fide assigned according to 
alignment consensus of the following panel. (B) Alignment of 
reference PR-2 sequences with JV. clevelandii peptides (shown above 
the alignment). Residues that align with the peptides identified are 
boxed; background according the ClustaLX color scheme [120]. In 
the rest of the alignment, only K/R and P are colored to 
highlight possible trypsin cleavage sites. Ntab_PR-2II, tobacco 
class II PR-2 isoform GI9 (GenBank P23547.1); Nben_PR-2II, 
JV. benihamiana class II PR-2 predicted protein (Sol Genomics 
NbC25742659g0003.1); Stub_PR-2II, Soknum tuberosum class II 
PR-2 (GenBank CAE52322.1); Slyc_PR-2II, S. lycopersicum class II 
PR-2 (GenBank NP_001234798.1); Ntab_PR-2I, tobacco class I 
PR-2 clone gglbSO (GenBank AAA63541.1) as oudier. The last 22 
residues of Ntab_PR-2I are hidden, as they do not align to the 
other accessions. In each alignment lane, positions of the first and 
last accession amino acids are indicated. In de novo defined 
peptides, residues that difiFer from the Ntab_PR-2II sequence are 
marked with a star. 
(TIF) 

Text SI This file contains: Supplemental Methods; 
Supplemental References; Table SI, Virus species and 
GenBank accession numbers for viral protein sequences 
used in the study; Table S2, List of the primers used in 
the study. 
(PDF) 
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